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1. INTRODUCTION 

Our theoretical understanding of superconductivity in magnesium diboride (MgB 2 ) has 
made rapid progress since its discovery by Nagamatsu et al. JH- Unlike superconductivity 
in the high-T c cuprates we are in a possession of a very clear picture of its superconducting 
state now \2\. It seems clear that superconductivity is driven by electron-phonon interaction in 
MgB 2 . More excitingly, it appears well established both theoretically and experimentally that 
the rare form of two gap superconductivity is realized in this compound. Two superconducting 
gaps of distinctively different size appear to exist on different disconnected parts of its Fermi 
surface. Since MgB 2 is the clearest example of two gap superconductivity to date, it makes 
it an interesting object for study and exploration. In the present work we want to review our 
present understanding of two gap superconductivity in MgB 2 from the theoretical perspective 
and discuss some of its consequences. The presence of these two different gaps gives rise to a 
number of anomalies and some of them shall be discussed here. Since the literature on super- 
conductivity in MgB 2 has grown rapidly, no complete coverage of experimental and theoretical 
work can be made here. Instead, the theoretical results most crucial to our understanding in the 
authors view will be discussed and a selection of peculiar consequences will be presented. 

Since much of our present understanding stems from band structure calculations and solu- 
tions of Migdal-Eliashberg equations, we are reviewing these results in the next section. The 
emergence and stability of two gap superconductivity in MgB 2 shall be discussed in section 3. 
In section 4 we will discuss consequences of this band structure picture for the upper critical 
field. Section 5 is finally devoted to the microwave conductivity. 



Figure 1 : The Fermi surface of MgB 2 . (Adapted from Ref. 0, ©2001 The American Physical 
Society). 



2. THE PICTURE SUGGESTED BY BAND STRUCTURE CALCULATIONS 

Magnesium diboride possesses a comparatively high critical transition temperature com- 
pared with other conventional superconductors of about T c =40 K. Presently only the high-T c 
cuprates have higher transition temperatures. For that reason the natural question arises whether 
superconductivity in MgB 2 is of the conventional electron-phonon driven type or if its super- 
conducting state has more similarities with the one in the high-T c cuprates. In the high-T c 
cuprates the superconducting state is of an unconventional d-wave type, possessing gap nodes, 
and electronic pairing mechanisms, like for example exchange of antiferromagnetic spin fluc- 
tuations, are discussed seriously. An ongoing discussion still concerns the role and importance 
of phonons in these compounds. 

In contrast to optimally doped high-T c cuprates MgB 2 shows a strong isotope effect. If the 
boron- 1 1 isotope is replaced by the lighter boron- 10 isotope, T c increases by about one Kelvin, 
indicating an important phonon contribution to the pairing interaction Q. Low temperature 
specific heat and penetration depth studies are consistent with an exponential decay, indicating 
the presence of a full gap without nodes BH|5]|. In addition, there are no indications of sizeable 
magnetic interactions in MgB 2 , again in contrast to the high-T c cuprates. All these observations 
seem to place MgB 2 among the conventional s-wave electron-phonon driven superconductors. 
Then the question arises, why it has such a high transition temperature as compared with other 
conventional superconductors. Here, band structure calculations turned out to be elusive, which 
we want to review in the following. 

The MgB 2 lattice structure consists of alternating layers of Boron and Magnesium atoms. 
The Boron atoms form a honeycomb lattice and the Magnesium atoms a triangular lattice 
halfway between the Boron layers. Calculations of the electronic band structure show four 



bands crossing the Fermi energy leading to four topologically disconnected Fermi surface sheets 
shown in Fig. d 1 6 1 . Two of these bands are derived from Boron p z orbitals. They form the so- 
called 7r bands seen as the red (electronlike) and blue (holelike) tubular networks in Fig.[T] The 
other two bands derive from Boron p^, and p y orbitals and form the so-called a bands, seen 
as the green and blue cylindrical Fermi surfaces centered around the T point (both holelike). 
These possess mainly two dimensional character. Interestingly, all these bands are dominated 
by Boron p orbitals and contributions from Magnesium orbitals are very small at the Fermi 
level. About 58 percent of the total density of states at the Fermi level is residing on the ir bands 
making both a and ir bands about equally important for the electronic properties of MgB 2 . 

Density-functional calculations of the phonon modes and the electron-phonon interaction 
strength can be found in Ref s . [ 7 1 , [ 8 1 , and [ 9 1 . The highest phonon density of states is found in 
the energy range around 30 meV. However, these phonons only couple weakly to the electrons 
at the Fermi level and thus do not contribute very much to superconductivity. This can be nicely 
recognized in Fig. 1 of Ref. Q, where the interaction strength of the phonons is shown as the 
area of the black circles in the figure. In fact, the phonons that couple most strongly to the 
electrons at the Fermi level are found in the energy range around 70 meV. These phonon modes 
evolve from the E 2g mode at the T point and correspond to a Boron-Boron bond- stretching 
vibration of the Boron sub lattice. A comparison with the phonon modes in the isostructural 
but nonsuperconducting compound A1B 2 in Ref. [8| shows that these E 2g phonon modes are 
strongly softened in MgB 2 consistent with their strong coupling. Correspondingly, the so-called 
Eliashberg function ol 2 F{uj), which weights the phonon density of states with the coupling 
strengths and appropriately describes the pairing interaction due to phonons, possesses a strong 
peak around 70 meV and significantly differs from the phonon density of states in contrast to 
conventional strong-coupling superconductors. The dimensionless electron-phonon coupling 
constant was found to lie between A ~ 0.7-0.9 from these first-principles calculations. From 
this microscopic information we can obtain a qualitative understanding of why T c is so high in 
MgB 2 by looking at the BCS T c formula (which is of course a bad approximation in the present 
case, but can give us some qualitative tendencies): 

k B T c = 1.13huj c e-^ VN ^ (1) 

Here, co c is a characteristic phonon frequency, V the interaction strength and iV(0) the density 
of states at the Fermi level, with A ~ VN(0). First of all, the characteristic phonon frequency 
is comparatively high, because Boron is a light element and the E 2g phonon modes in question 
only involve vibrations of the Boron sub lattice. Secondly, this high frequency phonon at the 
same time possesses a strong coupling to the electrons at the Fermi level. This means that in 
MgB 2 we have a favorable coincidence of two effects helping to raise T c . 

3. TWO GAP SUPERCONDUCTIVITY IN MAGNESIUM DIBORIDE 

It was noted already early on that there is a problem with the superconducting gap size in 
MgB 2 . From the BCS gap ratio A = 1.76A;bT c one should expect a gap of 6 meV or some- 
what larger, if strong electron-phonon coupling effects are considered. However, experimental 
results obtained from different experimental techniques seemed to scatter between 2 meV and 
8 meV with some clustering around 2.5 meV and 7 meV |[T0l . Initially, one might have at- 
tributed this to insufficient sample quality, however, by the end of 2001 high quality single 



crystals became available and this problem remained. Also, values of the coherence length £ 
extracted by different means turned out to differ considerably. For example, STM tunneling 
measurements of the vortex core size by Eskildsen et al. ifTTl were consistent with a value of 
£ ^50 nm, while the coherence length as determined from the upper critical field value on the 
same sample gave £ ~10 nm. Since the coherence length is related to the superconducting gap 
via £ = hvp /7rA these differences can be related to differences in the gap value as well. These 
problems suggested studying the possibility of anisotropic or multiple gap structures in MgB 2 . 

Again, important insight into this question was provided by Migdal-Eliashberg type cal- 
culations based on band structure calculations @ E21 EU. Migdal-Eliashberg theory is a 
generalization of BCS-theory and allows to take into account the detailed properties of phonons 
and their coupling to the electrons at the Fermi level in the pairing interaction. Decomposing 
the electron-phonon coupling constant into contributions from the four bands Liu et al. |9 1 were 
able to show that the pairing interaction differs considerably on the a and n bands. The pairing 
interaction on the a bands turns out to be much stronger than the one on the n bands and the 
interband pairing strength. The reason for this is that the E 2g phonon mode being an in-plane 
vibration of the Boron atoms preferentially couples to the in-plane electrons on the more two di- 
mensional a bands, while its coupling to the three dimensional 7r bands is much weaker. Due to 
this anisotropy the resulting effective coupling constant was shown to increase to about A ~ 1, 
being quantitatively more consistent with a T c of 40 K. 

A fully momentum dependent solution of Migdal-Eliashberg equations was presented by 
H. J. Choi et al in Refs. 1~i~3l and ifbH . In this calculation the fully momentum dependent band 
structure, phonon modes, and electron-phonon couplings obtained from first principles density 
functional calculations were used as a starting point to solve Migdal-Eliashberg equations in 
the full Brillouin zone. The Coulomb pseudopotential was taken as a constant /x* = 0.12, 
which reproduces the experimental T c . Neglecting impurity scattering, this procedure allowed 
to calculate the value of the superconducting gap at each point on the Fermi surface. The result 
is shown in Fig. O where the size of the gap is shown color coded on top of the Fermi surface 
structure. This picture impressively demonstrates the presence of a large gap value around 7 
meV on the a bands and a small gap value around 2 meV on the ix bands. 

Consistent with Ref. [9| the authors were able to demonstrate the important role of both 
anisotropy and anharmonicity of the phonons on T c . When anisotropy was neglected and only 
the isotropic Migdal-Eliashberg equations were solved, T c dropped from 39 K to 19 K. When 
anharmonicity on the phonon frequencies was neglected but anisotropy was kept, T c was seen 
to increase from 39 K to 55 K [13]. This shows that phonon anharmonicity is harmful for 
superconductivity in MgB 2 , but anisotropy is helpful, and for a quantitative understanding of 
superconductivity in MgB 2 both are needed. Choi et al. were also able to show that the isotope 
exponent observed experimentally could be quantitatively reproduced, if the anharmonicity of 
the phonons was taken into account. Also, the temperature dependence of the specific heat from 
these calculations was in good agreement with the experimental results. All these results seem 
to show that we can obtain a complete understanding in terms of two gap superconductivity in 
MgB 2 based on the picture shown in Fig. O 

An important question in this picture concerns the role of impurities and in particular inter- 
band impurity scattering, however. Conventionally one would expect that impurity scattering 
tends to equalize all gap values on the Fermi surfaces thereby reducing T c . This is one of the 
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Figure 2: The energy gap distribution on the Fermi surface of MgB 2 (color coded). (Adapted 
from Ref. OH, ©2002 Nature Publishing Group). 



reasons why two gap superconductivity or anisotropic s-wave superconductivity is rarely ob- 
served. This does not seem to be the case in MgB 2 . For example, samples with very different 
residual resistivities are observed to have essentially the same critical temperature T c . This is a 
behavior that one would expect in a single gap s-wave superconductor due to Anderson's theo- 
rem. However, in a two gap superconductor a strong change of T c with impurity concentration 
should be expected [fT51l. A very good answer to this question has been given by 1. 1. Mazin et 
al. in Ref. 1 16 1: due to the particular electronic structure of MgB 2 the electronic wave func- 
tions on 7r bands and a bands possess different parity symmetry. The n bands deriving from 
the Boron p z orbitals are mainly antisymmetric with respect to the Boron plane, while the o 
bands deriving from both Boron p x and p y orbitals are mainly symmetric. This disparity of a 
and 7r bands makes the impurity scattering matrix element between these two types of bands 
exceptionally small as compared with the impurity scattering matrix element within each of 
these bands. Using density functional supercell calculations for various impurities Mazin et al 
were able to show that the interband scattering rate appears to be one to two orders of mag- 
nitude smaller than the intraband scattering rates due to this band disparity. This means that 



impurity scattering will equalize the gaps within each of the two types of bands, but not very 
much between them. This makes the two gaps in MgB 2 exceptionally stable against impurity 
scattering. Only a very large amount of impurity scattering is expected to lead to a reduction of 
T c due to interband scattering and an accompanied merging of the two gaps. 

While signatures of the two gaps have been observed in several different experiments, so 
far direct experimental confirmation of the important role of the E 2g phonon for the pairing is 
still lacking. In principle, such a strong coupling phonon could produce observable structures 
in the tunneling density of states. However, there are difficulties in the case of MgB 2 as has 
been discussed by Dolgov et al ifPTl : on one hand the high frequency of the phonon mode leads 
to a strong reduction of the structures produced in the tunneling spectrum, on the other hand 
tunneling spectroscopy is mostly dominated by the n band phonons, which are not so important 
for the pairing interaction. This makes observation of the E 2g phonon and its coupling strength 
in tunneling spectroscopy difficult. 

Concluding this section we can say that we seem to have a consistent and quantitative pic- 
ture of two gap superconductivity in MgB 2 , which mainly arises from density functional cal- 
culations and solutions of Migdal-Eliashberg equations. The two gaps arise due to a strongly 
anisotropic electron-phonon interaction of the (anharmonic) E 2ff phonon mode, which couples 
more strongly to the electrons on the a bands than to the ones on the n bands. Due to the dispar- 
ity of the two types of bands these two gaps appear to be very stable against impurity scattering, 
allowing T c to remain large even in samples with large residual resistivities. 

4. UPPER CRITICAL FIELD ANISOTROPY 

In this section we want to explore the consequences of the two gap picture presented in 
the previous section on the temperature dependence of the upper critical field B c2 , particularly 
its anisotropy. Measurements of the upper critical field in MgB 2 single crystals have shown 
that not only B c2 is anisotropic, but also that this anisotropy is strongly temperature dependent 
lfT8lfT9ll20ll2Tl . The anisotropy ratio r c2 is given by T c2 = B^ 2 /B^ 2 , where B^ 2 is the upper 
critical field, when the field is applied along the Boron planes, and B c c2 the one for field along 
c-axis direction perpendicular to the Boron planes. At low temperatures T c2 reaches values 
around 5. It decreases with increasing temperature and reaches values around 2 at T c (see the 
solid circles in Fig. |4] below). Close to T c B c c2 varies linearly with temperature, while B% 2 
shows a pronounced upward curvature. This behavior is quite unusual, because in conventional 
single gap superconductors the anisotropy ratio is very much temperature independent, rarely 
changing by more than 10 to 20 percent. An upward curvature of the upper critical field has 
been observed in YNi 2 B 2 C and attributed to two band behavior, however, in this case the upward 
curvature appeared in all spatial directions, in contrast to MgB 2 ll22l . 

The first theoretical work addressing this unusual behavior in MgB 2 was a study within an 
anisotropic gap model by Posazhennikova et al. In this work a single anisotropic s-wave 
gap on an elliptical Fermi surface was considered. It was shown that the experimental data of the 
upper critical field in MgB 2 single crystals including both the strong temperature dependence 
of the anisotropy ratio as well as the upward curvature only appearing with field in afe-plane 
direction could be fitted within this model, if the gap was smaller in c-axis direction than in ab- 
plane direction. While this model can account for all experimental observations on the upper 
critical field in MgB 2 , it possesses some drawbacks, however, regarding its consistency with 



other experimental results. One problem concerns the ratio of the maximum to the minimum gap 
value. Within this anisotropic s-wave model a gap ratio of about 1 : 10 was needed to explain the 
upper critical field data. This ratio seems too large compared with the about 1:3 ratio as observed 
in tunneling experiments, for example. In addition, measurements of the in-plane penetration 
depth clearly show an exponential decrease with the small gap ll24ll25ll . This suggests that the 
small gap is present within the Boron plane direction consistent with the picture in Fig. |3 but 
inconsistent with this anisotropic s-wave model. 

Another problem concerns the anisotropy of the lower critical field B c \. From anisotropic 
(single gap) Ginzburg-Landau theory one should expect that the upper critical field anisotropy 
is related to the lower critical field anisotropy via B^jB c c2 = B^/B^. Recent experimen- 
tal studies have established, however, that this relation is violated in MgB 2 Il26ll27ll28 l. the 
upper critical field anisotropy decreasing with increasing temperature, while the lower critical 
anisotropy is found to increase. This is in agreement with expectations based on two gap models 
1 29], but again inconsistent with an anisotropic s-wave model. 

A first calculation within an anisotropic two band model was presented by Miranovic et 
al. ll29l based on clean limit Eilenberger theory. The Fermi surface sheets of the two bands 
were taken as ellipsoids with very different anisotropics, while the two gaps were taken to 
be isotropic. Within this calculation it was also possible to qualitatively reproduce the ex- 
perimentally observed upper critical field anisotropy as in the anisotropic s-wave model by 
Posazhennikova et al. In addition to that, however, this two band calculation did not suffer from 
the drawbacks mentioned above: more realistic values for the ratio of the two gaps could be 
taken, the small gap was present within the Boron plane direction, and the lower critical field 
anisotropy was found to increase with temperature. 

In order to answer the question whether the picture in Figgis consistent with the upper 
critical field anisotropy the author and N. Schopohl presented a calculation within two band 
Eilenberger theory taking into account the detailed band structure Fermi surface topology based 
on Fig. [T] [ 30]. For that purpose a simple but realistic model of the Fermi surface was used as 
shown in Fig. |3] The upper panel of Fig. |3] shows one half of the Brillouin zone from Fig. [T] 
Using the periodicity of the lattice identical copies of the Brillouin zone may be attached to each 
other (center panel). Choosing a new Brillouin zone around the T point leads to the structure 
shown in the lower panel. This rearrangement shows that the a band Fermi surface can be 
modeled by a distorted cylinder and the n band Fermi surface by a torus cut open half. Here, 
the two a bands and the two n bands are assumed to be degenerate, for simplicity. This suggests 
the following parameterization of the Fermi surfaces for the a band Fermi momentum OTI: 
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and the n band Fermi momentum: 
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Figure 3: Simplification of the MgB 2 Fermi surface from Fig.d Rearranging the Fermi surface 
sheets around the T point shows that the a band can be modeled as a distorted cylinder and the 
7r band as a half-torus. 



Here, <p is the in-plane angle, k z the c-axis momentum, and d the azimuthal angle of the half- 
torus running from 7r/2 to 3tt/2. The corresponding Fermi velocities are given by 
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Here, the dimensionless c-axis dispersion parameter e c has been assumed small. For MgB 2 the 
following parameter values have been extracted from band structure calculations ll30ll : v FjCr = 
4.4 x 10 5 m/s, v F)V = 8.2 x 10 5 m/s, e c = 0.23, and k = 0.25. 

Based on this model for the MgB 2 Fermi surface one can now calculate the upper critical 
field using the linearized two band Eilenberger equations. For u n > they read 



jcL> n + V F ,a 

along with the gap equation 

A a (r) = 
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Here, f a is the anomalous Eilenberger propagator and A a (r) the (spatially dependent) gap 
function for the two bands a E {n,a}. The pairing interaction \ aa becomes a two-by-two 
matrix in the band indices. The effective matrix elements have been calculated from band 
structure calculations in Ref . [ 9 1 and found to be 
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Together Equations © and © establish an eigenvalue problem for the gap function A Q (r). For 
a given temperature T one has too look for the solution A a (r) that solves both equations for the 
highest value of the magnetic field B = V x A. 

Eilenberger theory is a generalization of BCS theory to inhomogeneous superconducting 
states and contains Ginzburg-Landau theory as a limiting case for T — > T c ll36l 1371 . It holds 
in the limit k F £ 1, where £ is the coherence length and k F the Fermi momentum. Since 
Ginzburg-Landau theory is limited to the vicinity of T c , Eilenberger theory is the method of 
choice, if one wants to calculate properties of type-II superconductors in the vortex state at 
lower temperatures from microscopic grounds. Near the upper critical field the gap function 
becomes small and therefore one can use the linearized Eilenberger equations to determine B c2 . 

Equations © and © are usually solved using a Landau level expansion of A Q (r) above 
the Abrikosov ground state of the vortex lattice [32]. For strongly anisotropic systems this 
procedure requires a high number of excited states, however. In addition, the calculation can 
only be done numerically. In Ref. [ 30 ] Equations © and © were solved using a variational 
method that had been introduced earlier by the author in Ref. [23|. This method does not 
require a high number of excited states and even allows to obtain analytical results in some 
limiting cases. The variational ansatz for A a (r) corresponds to a distorted Abrikosov ground 
state of the form 

(8) 
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Figure 4: Temperature dependence of the anisotropy ratio T c2 

based two-band model of Ref. QUI . Results are shown for different interband pairing strengths 
rj and c-axis dispersions e c . Solid circles are experimental results from Lyard et al. [20] (adapted 
from Ref. EDI, ©2003 The American Physical Society). 



where r is used as a variational parameter and determined such as to maximize the upper critical 
field B c2 . Here, ip^ is the usual (undistorted) Abrikosov ground state. More details about this 
variational method can be found in the appendix. 

Using this method the temperature dependence of the anisotropy ratio was calculated from 
microscopic grounds based on band structure calculations in Ref. [ 30 ] as pointed out above. 
The result is shown in Fig.|U Here, the parameter r] is a dimensionless parameter describing the 
interband pairing strength and is given by 

\ 7T7V \ 

where A + and A_ are the larger and smaller eigenvalues of the matrix X aa ' in Eq. ©, respec- 
tively. The dashed line shows the result for the parameters obtained from band structure calcu- 
lations, the solid circles are experimental results on MgB 2 single crystals from Lyard et al. ||2"01 . 
For comparison, also the result with no interband pairing 77 = is shown as the solid line. In 
this case the upper critical field anisotropy is determined by the a band only. 

In Ref. O0l it was found that the two most important parameters determining the tempera- 
ture dependence of the upper critical field anisotropy are r\ and the c-axis dispersion e c of the 
cylindrical a band. All other parameters have only minor influence on the shape of the curves 
in Fig. |4j If these two parameters are allowed to vary from their band structure values, one 
can obtain an excellent fit of the experimental data shown as the dotted line in Fig. |4] In this 
(two parameter) fit the interband pairing strength 77 turns out to be somewhat larger and the 
c-axis dispersion e c a little bit smaller than expected from band structure calculations. Simi- 
lar observations, that larger interband pairing strengths than found by band structure calcula- 
tions are needed to fit experimental data, have been made also in Refs. tl33ll and ll34l. Possible 
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Figure 5: Anisotropy ratio T c2 = B^jB c c2 as a function of interband pairing strength i] for the 
band structure based two-band model of Ref. OOll at T = (dashed line) and T = T c (solid 
line) (adapted from Ref. |30|, ©2003 The American Physical Society). 



explanations could be either that band structure calculations underestimate the interband pair- 
ing strength or that some small amount of interband impurity scattering or effects of strong 
electron-phonon coupling, which were neglected in the above calculations of the upper critical 
field, mimic a somewhat stronger interband pairing strength. 

The result in Fig. |4] shows that indeed a strong temperature dependence of the upper critical 
field ratio consistent with experimental results has to be expected in the clean limit for the 
parameters found in band structure calculations for MgB 2 . As the solid line shows, the influence 
of both gaps is crucial here. The strong temperature dependence can be understood as follows: 
at low temperatures and high magnetic fields the a band with the large gap is dominating leading 
to a strong anisotropy. When temperature is increased and thus the magnetic field reduced, the 
influence of the n band with the small gap becomes more and more important. Since the n band 
is more isotropic a reduction of the upper critical field anisotropy results. 

It is instructive to look at the change of T c2 as a function of the interband pairing strength 
r] at T = and T = T c . This is shown in Fig. When r] = there is no strong temperature 
dependence of T c2 . When r\ is increased, T c2 decreases much more rapidly at T = T c than 
at T = 0, because the influence of the n band is more important at higher temperature. As a 
result the temperature dependence becomes stronger. However, when r\ is increased towards 0.5 
(maximum coupling of the two bands), eventually Y c2 becomes small also at low temperature 
and the temperature dependence becomes weak again. It appears that either for very weak or 
for very strong interband pairing interaction the system effectively behaves like a single gap 
superconductor with a weak temperature dependence of the upper critical field anisotropy. It is 
only in the intermediate regime around rj ~ 0.05 — 0.3 where a strong temperature dependence 
of T c2 can be observed. This is just the regime in which MgB 2 appears to be according to band 
structure calculations. 

In Ref. 031 further results of this model can be found: the temperature dependence of the 



upper critical field for field applied along the Boron plane direction shows an upward curvature, 
while for field applied in c-axis direction it does not, in good agreement with the experiments. 
Also, the magnetic field dependence of the zero energy density of states follows the experimen- 
tal temperature dependence of the electronic specific heat coefficient. 

The theories about the upper critical field in MgB 2 discussed so far are clean limit theories, 
e.g. impurity scattering was assumed to be small. However, even in the best MgB 2 single crys- 
tals available so far it is believed that in the n band the scattering rate is larger than the gap, 
while in the a band the scattering rate might be smaller than the gap. This means that impurities 
are expected to affect the temperature dependence of the upper critical field and one might ask 
how this could affect the clean limit results discussed above. So far no calculations for general 
impurity scattering rates have been presented. However, calculations can be considerably sim- 
plified in the dirty limit. In Refs. Il3lfll3*9l such dirty limit theories of the upper critical field for 
MgB 2 were presented. The main simplification here is that the intraband scattering rates rV and 
T a are assumed to be larger than the gaps, while the interband scattering rate is assumed to be 
negligible due to the parity argument by Mazin et al. [ 16 1 discussed in the previous section. In- 
terestingly, the strong temperature dependence of the upper critical field anisotropy was shown 
to exist also in this dirty limit [39|. The physical reason for this behavior is the same as in the 
clean limit: an interband pairing interaction of intermediate strength and a strong anisotropy of 
the Fermi velocities in the a band which results in a strong anisotropy of the diffusivities in this 
band. Therefore, the inclusion of impurities qualitatively does not change this behavior. 

In Ref. [39 1 also the angular dependence of the upper critical field was studied. It was 
shown that this angular dependence shows strong deviations from the results expected from 
anisotropic Ginzburg-Landau theory. The reason for this is that the c-axis coherence lengths in 
the two bands strongly differ. As a result the validity of Ginzburg-Landau theory is reduced to 
a very narrow region near the critical temperature T c iRUl . 

To summarize this section we have seen that the unusually strong temperature dependence 
of the upper critical field anisotropy can be nicely understood in terms of the two gap model, if 
the particular Fermi surface topology of MgB 2 is taken into account. The essential ingredients 
for this effect are two Fermi surfaces with very different anisotropics and an interband pairing 
interaction with an intermediate strength. Too weak or too strong interband pairing would result 
in effectively single gap behavior. According to band structure calculations MgB 2 fulfils these 
requirements possessing a strongly anisotropic cylindrical o band and a more isotropic toroidal 
7r band. Parameters from band structure calculations can even reproduce the upper critical field 
anisotropy quantitatively. This shows that the behavior of the upper critical field is consistent 
with the band structure picture presented in the previous sections. 

5. MICROWAVE CONDUCTIVITY 

In this section we want to discuss the consequences of two gap superconductivity on the 
microwave conductivity in MgB 2 , particularly its temperature dependence at fixed frequency. 
In Ref. [33 1 the observation of an anomalous microwave conductivity peak in MgB 2 thin films 
has been reported and its interpretation in terms of two gap behavior shall be discussed here. 
We will first start with a review of conductivity peaks in conventional superconductors as well 
as in high-T c cuprates and then compare them with the results in MgB 2 . 

In conventional dirty limit superconductors the microwave conductivity at frequencies suf- 
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ficiently below the gap frequency shows a coherence peak as a function of frequency, which 
is related to the so-called Hebel-Slichter peak in the temperature dependence of the NMR spin 
relaxation rate. As shown in Fig.|6]the conductivity initially increases when entering the super- 
conducting state going through a peak value near 0.9 T c and finally being suppressed exponen- 
tially at lower temperatures. The peak naturally comes out of BCS theory and at the time was 
regarded as one of the key triumphs of BCS theory, because earlier theories based on two fluid 
models were not able to account for this effect. When the superconducting state is entered and 
the gap opens, a square-root singularity appears in the density of states at the gap edge. This 
singularity leads to increased contributions to the microwave conductivity. When the gap A 
increases upon lowering the temperature, it eventually becomes larger than k B T. At this point 
quasiparticles are condensed out into the superfluid and the microwave conductivity is sup- 
pressed exponentially at lower temperatures. As a rule of thumb we can estimate the position 
of the coherence peak roughly at the point where A(T) becomes equal to k B T, which is shown 
as the crossing point between the solid and the dashed line in the inset of Fig. |6] In contrast 
to the NMR Hebel-Slichter peak, the coherence peak in the microwave conductivity depends 
on the scattering rate, however. When the system becomes cleaner, the microwave coherence 
peak is gradually reduced, as has been discussed by Marsiglio [41 1. The ultimate reason for this 
dependence on the scattering rate is lying in the fact that the conductivity is a nonlocal quantity 
in contrast to the NMR relaxation rate which is a local probe. 

In the high-T c cuprates the situation is completely different. In measurements of the NMR 
relaxation rate no Hebel-Slichter peak is observed. However, the microwave conductivity shows 
a huge peak at temperatures between 0.3 T c and 0.6 T c depending on frequency as shown in 
Fig. [7] The absence of the Hebel-Slichter peak can be easily understood as a consequence of 
the (i-wave nature of the superconducting state. In a ci-wave superconductor the singularity at 
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Figure 7: Temperature dependence of the microwave conductivity in YBCO showing a large 
conductivity peak (adapted from Ref. ©1999 The American Physical Society). 



the gap edge of the density of states is not a square-root singularity anymore, but a logarithmic 
singularity, having a much weaker influence on the NMR relaxation rate, as shown in Fig. [8] 
In addition, strong-coupling effects tend to wash out the singularities in the density of states, 
which leads to an additional suppression of the coherence peak as discussed in Ref. iflTl . 

Since there is no coherence peak in the NMR relaxation rate, why is there a peak in the 
microwave conductivity? It has been suggested that this peak has a different physical origin and 
is related to a rapid suppression of the inelastic scattering rate in the superconducting state [1431. 
Such a drop of the inelastic scattering rate naturally appears when inelastic scattering in the 
normal state is dominated by electron-electron scattering, for example by spin fluctuations, be- 
cause the superconducting gap is suppressing electron-electron scattering below the gap energy 
in the superconducting state. Detailed theoretical calculations by Hirschfeld et al l44ll including 
(i-wave superconductivity, strong impurity scattering, and inelastic spin fluctuation scattering 
were shown to be consistent with the experimental results. 

In MgB 2 microwave conductivity measurements have shown a conductivity peak appear- 
ing at about 0.5 T c as shown in Fig.0 This peak position seems to be somewhat intermediate 
between conventional superconductors and high-T c cuprates and the natural question arises, 
whether this is a coherence peak or a peak due to lifetime effects as in the cuprates. In order 
to address the question of lifetime effects a quick look at the temperature dependence of the 
resistivity in MgB 2 is useful. It has been shown that the temperature dependence of the resis- 
tivity can be fitted well by Bloch-Griineisen formula with a temperature independent residual 
resistivity associated with impurity scattering and a temperature dependent phonon scattering 
part [45|. This seems to show that no large electron-electron scattering is present in the normal 
state of MgB 2 . Moreover, at T c when superconductivity sets in, the resistivity is already in the 
saturation regime where phonon scattering is frozen out (due to the high frequency of the E 2g 




Figure 8: (a) Density of states in a rf-wave superconductor (solid line) and an s-wave super- 
conductor (dashed line), (b) Corresponding NMR relaxation rates for <i-wave (solid line) and 
s-wave (dashed line). In the <i-wave case no large Hebel-Slichter peak is visible anymore. 




Figure 9: Temperature dependence of the microwave conductivity in MgB 2 thin films (cir- 
cles and squares). For comparison, the microwave conductivity of a Niobium thin film in the 
dirty limit is also shown (triangles). (Adapted from Ref. P51 . ©2003 The American Physical 
Society) 



(a) (b) (c) 

Figure 10: Diagrammatic contributions to the conductivity, (a) The main contribution to the 
7r band, (b) The main contribution to the a band, (c) The lowest order interband scattering 
contribution to the conductivity which is expected to be small in MgB 2 . 



phonon mode). Therefore the resistivity at T c seems to be dominated by elastic impurity scatter- 
ing and thus no rapid drop of inelastic scattering can be expected in the superconducting state. 
Note that this is in contrast to the high-T c cuprates, where resistivity varies linearly with tem- 
perature down to T c and is not saturated. These arguments seem to exclude the interpretation of 
the conductivity peak in MgB 2 in terms of a lifetime effect. 

In the following we therefore want to study the microwave conductivity of a dirty two gap 
superconductor. Dirty in the sense that the intraband impurity scattering rates in the two bands 
are assumed to be larger than the respective gaps, but interband scattering is neglected as sug- 
gested by the argument by Mazin et al. ifTol . 

In a two band system like MgB 2 we can approximate the total conductivity as the sum of the 
two partial conductivities of the two bands (parallel resistor formula). This approximation ne- 
glects interband scattering events like the one shown diagrammatically in Fig. fToT c). However, 
these events are expected to give only small corrections, because interband impurity scattering 
is expected to be much weaker than intraband impurity scattering due to the argument by Mazin 

etal. ma. 

If the intraband scattering rates in the two bands are sufficiently larger than their gaps, we 
can use the Mattis-Bardeen dirty limit formula for the conductivity in each band separately 

aJu) 1 r°° ( , Q + uj , Vt 



roc / () _|_ rj () \ 

/ dQ tanh — tanh — [N a (fl) N a (ft + u) + M a (fi) M a (Q + u)} 

J — oo 

(10) 



&n,a 2^; J- 

Here, the normal and anomalous densities of states are given by 




N a (tt) = Re{ j=p = } and M a (Q) = Re { A "^ gn B } (11) 



The total conductivity under these assumptions is given by ct\{uj) = cr^iu) + a a (uj). The 
temperature dependence of the two gaps A Q has to be determined from a solution of the two by 



0.5 1 

t=TA c 



Figure 1 1 : Temperature dependence of the two gaps in MgB 2 found from a solution of the two 
by two gap equation. The dashed line shows A = k B T. Its crossing point with the small gap 
appears at a much lower temperature than the one with the large gap. 



two gap equation: 



Here, uo c is a characteristic phonon cut-off frequency and \ a " the two by two coupling matrix. 
These parameters can be either taken from a band structure calculation or tried to be adjusted to 
the particular sample in question. In Ref. an intermediate approach was taken: the partial 
densities of states and the cut-off frequency were taken from band structure calculations, while 
the remaining parameters were adjusted to the T c of the sample and the value of the small gap, 
which could be extracted from the exponential fall-off of the temperature dependent change 
of the penetration depth. A typical temperature dependence of the two gaps found this way is 
shown in Fig.[TT] One notices that the small gap reaches k B T at a much lower temperature than 
the large gap. According to the rule of thumb given above this would mean that the coherence 
peak in the tc band is expected to appear at much lower temperature than in a conventional 
superconductor. 

That this is indeed the case is shown in Fig. [121 Here, the temperature dependence of the 
Mattis-Bardeen conductivity is shown for different zero temperature gap values. In these calcu- 
lations the BCS temperature dependence has been scaled to different gap values for simplicity. 
Clearly, the position of the coherence peak follows the crossing point of the gap with k B T. 

In Fig.|9]the fits to the two MgB 2 samples show the temperature dependence of the n band 
contribution to the conductivity obtained using the temperature dependence of the small gap in 
Fig. CU Apparently there is good agreement between measurement and calculation. The fit to 
the Niobium sample was obtained using the BCS temperature dependence of the gap, as one 
should expect for a conventional superconductor. 

Now one might ask where the contribution from the a band in the MgB 2 samples is. Accord- 
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Figure 12: Temperature dependence of the microwave conductivity (upper panel) for different 
gap values (lower panel) at a frequency of uo = 0.02k B T c . 
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Figure 13: Temperature dependence of the microwave conductivity for different ratios of the 
scattering rates r <T /r 7r =0.1, 1, and 10 along with the experimental results from Ref. Q3l. 



ing to Fig.[l2]the o band should produce a conventional coherence peak. The relative weight of 
the two contributions is determined by the partial normal state conductivities of the two bands, 
however. We have 



= a w (u) + cv(cj) = a n ^— h o n ^— (13) 

where the conductivity ratios are given by Eq. dTOb . The normal state conductivities are related 
to the plasma frequencies u) p>a and the intraband scattering rates r a via the expression 

u 2 

<?n,a = ^O^T (14) 

Theoretical band structure values for the plasma frequencies have been given in Ref. IIT6 I. Ac- 
cording to these calculations we have cu Ptn = 5.89 eV and uj p ^ = 4.14 eV for the in plane 
plasma frequencies. Therefore, if the scattering rates in the two bands are the same, the con- 
tribution of the 7r band is expected to be a factor of 2 larger than the one of the a band. This 
is not sufficient to understand the apparent absence of a sizeable a band contribution in the 
experimental data. A possible explanation could be a stronger scattering rate in the a band at 
least in the thin film samples studied in Ref. ll33ll . To illustrate that, in Fig. [13] the temperature 
dependence of the conductivity is shown for different values of the relative ratio of the two scat- 
tering rates. These results suggest a at least 10 times larger scattering rate in the a band than in 
the 7r band. In principle, the scattering rates in the two bands can be varied by selective doping, 
for example by doping with Aluminum for Magnesium or Carbon for Boron ll48ll49lfT6ll . It 
would be interesting to see whether a cross-over like the one shown in Fig.[l3]can be observed 
experimentally on a series of MgB 2 films with varying degree of doping. 

Concluding this section we have seen that the anomalous microwave conductivity peak ob- 
served in MgB 2 thin films can be understood as a coherence peak due to the small gap. The 
peak appears at much lower temperatures than in conventional superconductors, because the 
small gap opens up more slowly upon reducing the temperature below T c . As a result the con- 
densation of excited quasiparticles and with it the exponential suppression of the conductivity 
sets in later. Again, this peculiar effect is consistent with two gap behavior and reinforces the 
picture presented in the previous sections. 

6. SUMMARY 

We have reviewed our present picture of superconductivity in magnesium diboride. Much 
of this picture stems from band structure calculations. We have seen that the high value of 
the critical temperature T c is due to a high frequency phonon mode that couples strongly to 
the electrons at the Fermi level. This coupling is particularly strong to the electrons on the o 
band and leads to a much larger superconducting gap on the a band than on the n band. The 
different parity of the electronic wavefunctions of the a band and the n band results in a strong 
reduction of the interband impurity scattering matrix element. This makes the two different gaps 
particularly stable against impurity scattering. For these reasons magnesium diboride appears 
to be the clearest example of an intrinsic two gap superconductor to date. 

Two gap superconductivity in magnesium diboride leads to unusual behavior in a number 
of experimentally accessible quantities. Here, we have focused on the temperature dependence 



of the upper critical field anisotropy and the microwave conductivity. The upper critical field 
anisotropy shows an unusual strong temperature dependence. This feature is related to the 
presence of the two gaps, but also to the very different anisotropy of the two bands. Here, 
the intermediate strength of the interband pairing interaction plays a crucial role as well. The 
temperature dependence of the microwave conductivity shows a peak at temperatures around 
0.5T C . This peak can be understood as an anomalous coherence peak that is shifted downwards 
in temperatures because of the small gap. 

In conclusion, the case of two gap superconductivity in magnesium diboride generated by 
strong electron-phonon coupling appears to be well established and can be consistently ob- 
served in very different experimental quantities such as the upper critical field or the microwave 
conductivity. 

The author would like to thank O. V. Dolgov, S. Graser, C. Iniotakis, B. B. Jin, N. Klein, 
S.-I. Lee, K. Maki, A. I. Posazhennikova, and N. Schopohl for valuable discussions about this 
and related topics. 



APPENDIX 



In this appendix the variational method for the calculation of the upper critical field from 
Eqs. © and © shall be described. This method was used in Refs. Il23ll30l . 
Defining the operator 



L a = 2 \u n \ + sgn u n v F:a (Pj 



' - 2e - 
hV -i—A(r) 
c 



(15) 



where VF,a(k) is the Fermi velocity of band a, A the vector potential due to the internal magnetic 
field within the system, and u n = (2n+l)nT the Matsubara frequencies, Eq. © can be inverted 
using the identity 

L~ l = / ds e~ sLa (16) 



which leads to 

poo 

f a (r, k; u n ) = -2 / ds e~ sL « A Q {f) (17) 



Introducing this into the gap equation Eq. © we can eliminate f a : 

A a (f)=27rT]T ]T A m '(f^e- si »'A a) (r)) (18) 

This equation is a linear equation for A a (f). 

As is usual in weak-coupling theory we can eliminate the cutoff frequency u c in favor of the 
critical temperature T c . For this purpose we consider Eq. (fT8l in the absence of a magnetic field 
at T c . Then the gap function A a {f) becomes homogeneous and we find 

roo 

A a = 2vrT c ^ / rfse- 2s|w " (Tc)l 53A Qa A a , 

|u>n(T c )|<w c a' 

= 2 « T c E -4rE Am ' A «' d9) 

u)n(T c )>0 Un ^ ±c > a' 

This is an eigenvalue equation for A a and the largest eigenvalue A + of the matrix X aa ' deter- 
mines T c . Thus we find 

1 fc 1 

lT = 2?rTc E TTfTT (20) 

For tu c 3> T c the following relation holds: 

2^c E — T^tt - 27rT £ — — -In- (21) 
and thus we can write Eq. d20b in the form 

J_ _i n ^ = 2 7rT £ _L_^ = 2vrT £ / dse"^^ (22) 

A + Tc ^(T)>0^( T ) [w n (T)|<w c 



Multiplying this equation by J2 a > A aQ 'A a /(r) it can be subtracted from Eq. (TT81) leading to 

Aa(r=) + E A° Q ' A a / (r) (-±- + In ^] 



a' |w n |<w c 
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(23) 



Here, we have eliminated the sgn u n factor assuming inversion symmetry of the Fermi velocity 



(k) 



~V F ,a ~ 



(—k) and introduced the operator 



i c 



(24) 



In Eq. d2"3l we may now extend the uj n summation to infinity and sum it up: 



E « 
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-1su n 
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-2s7rT 
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1 - e- 4swT 2 sinh (2svrT) 



(25) 



Note that integration and summation in Eq. d23t may only be interchanged because the diver- 
gence for s — > has been eliminated. (In Eq. (TT81 this summation was not possible). So we 
find from Eq. 



Aa(f) + E A aQ 'A a ,(r) + In 

a' \ + 



T 



o sinh u 



e -iuv Fa ,-n/(27rT) _ 1 
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(26) 



In the presence of an external magnetic field B = V x A it is convenient to choose the field 
direction as the z-axis of the coordinate system. In these coordinates A a (r) does not depend on 
z. Choosing the gauge A = Bxy the operator n simplifies to 



-> h-i 2e i , . 
n = -V A (f) 

i c 



V 



-ihd x 




-ihd y - ^Bx | = \l — 



leB 
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( a + a) 



t{a — a 




(27) 



where we have introduced raising and lowering operators 2y ^-a = —ihd x — hd y + 2i^-x and 

^\/? a ^ = — ifcdx + hdy — 2i^-x. Thus, we see from Eq. d2"6l) that only the components of the 
Fermi velocity perpendicular to the magnetic field direction play a role and we have 



v FyCt • n = \ — 



K v x>ot + iv y>a ) a + (v x>a - ivy^ a ) a T 



(28) 



where v x a and the components of the Fermi velocity perpendicular to the magnetic field 

B. Note, that these components are not necessarily identical to the components of the Fermi 
velocity with respect to the crystal axes in Eq. ©. 

It is useful to introduce a scaling of the x- and y-coordinates of the form 



x = e x y — e y 



(29) 



Here e T is a scaling factor which scales the x- and y-coordinates in such a way as to preserve 
the area. Expressed in the new coordinates x and y the operator II can be written: 



n 



— ihe T <9 5 




-ihe T dy - ^e T Bx 




(30) 



where the raising and lowering operators a and are also expressed in terms of these new 
coordinates. Using this result we find 



v F>a • II = W — [(e v Xja + 



ic Vy^ a I a + | e v a 
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We can now write the exponential operator in Eq. (l26l) as follows: 



e -iuv F , a -Tl/(2irT) _ e A+B 



with 

A = 
B = 

Using the identity 

and calculating the commutator 
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we can write 
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(33) 
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(35) 



In principle, we can view Eq. (|2"6T) as an eigenvalue problem. The highest eigenvalue and 
its corresponding eigenfunction of the right hand side operator will determine the upper critical 
field B c2 (T). We already know from Abrikosov's solution of Ginzburg-Landau theory that the 



solution for an isotropic s-wave superconductor will be Abrikosov's vortex state wave function 
il)h{r). This wavefunction has the property that it is destroyed by the operator a 



a ^ A (f) = 



(36) 



and has the form of the lowest Landau level wavefunction. In principle, one can try to solve 
Eq. (l26t by making a Landau level expansion of A a (r) [32 , 50|. However, for an anisotropic 
superconductor one should expect a distortion of the vortex lattice. Therefore here we will start 
from a variational ansatz to Eq. (l26t by choosing a different wavefunction, which obeys the 
equation 

a A a (f) = (37) 

This corresponds to the choice Eq. (HI). Here, we can use the scaling parameter r now as a 
variational parameter which has to be adjusted such as to yield the highest possible eigenvalue 
of Eq. (l26b . If we insert this variational wavefunction into Eq. (l26t the e B and e A operators just 
drop out and we are left with the equation 
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where the function l a is given by the expression 
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(39) 



Equation d3*Ht is a 2 x 2 matrix equation and the upper critical field is determined by the cri- 
terion that the largest eigenvalue of the right hand side becomes 1. This criterion leads to the 
characteristic equation 



(1 - rj)l a + 7 1 l 7T + \nt = -A± (l a + hit) (Z ff + lnt) 
where t = T/T c and the parameters r\ and A± are given by 

and A ± 
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(40) 
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Here, A + and A_ are the larger and smaller eigenvalue of the matrix X aa ' , respectively. Note, 
that Eq. (ROT ) means that only the two parameters i] and A± out of the four parameters \ aa ' 
determine the upper critical field of a two gap superconductor. 

Equation (l40l) is a quadratic equation in hit and allows to calculate t once l w and l a are 
known. Therefore, calculation of B c2 from Eqs. (l39l) and d40l can proceed as follows: for given 
values of B c2 /T 2 and r the quantities l n and l a are calculated using Eq. (l39l) and the Fermi 
velocities given in Eq. ©. Then t can be calculated using Eq. (l40l) and the matrix elements 
given in Eq. ©. Now, the parameter r is optimized such as to maximize t. From the maximized 
t and the given value of £> c2 /T 2 finally B c2 can be calculated. This procedure is repeated for 
several values of B c2 /T 2 producing a B c2 {T) curve. 



In Eq. (l39t the integration over u has a special form and it is useful to introduce the integral 

r°° du 
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(42) 



While this integral cannot be calculated analytically in the most general case, it at least possesses 
some simple limiting expressions for small and large arguments y. (Large y corresponds to the 
zero temperature limit in Eq. ( 13*91 ). while small y corresponds to T — > T c , when B c2 — > 0.) 
These limiting expressions read 



K(%-fC(% 2 for y«l 



\ ln(4 7 y) 



for )/ > 1 



(43) 



Here, is the Riemann Zeta function and In 7 = 0.577215 Euler's constant. Using the 
integral I(y) Eq. (l39l) can be written 
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Using Eq. d43t this allows to find simple limiting expressions for l a in the limits T — ► and 
T — > T c . In particular, in the limit T — > T c we have in linear order in B c2 
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(45) 



Since B c2 goes to zero for T — > T c , both l a and In t decrease linearly and the right hand side of 
Eq. ("""01 can be neglected, because it becomes quadratic in 1 — t. Thus in linear order we have 



>4 



7g3)eB c2 

16C7T 2 T 2 



• {e~ 2T [(1 - V) + V «.)J + [(1 - V) + V «„) J } (46) 



Minimizing this expression with respect to r one finds 



-,2r 



Using this result we can calculate the slope of B c2 at T c from Eq. (|4"""1) : 



(47) 
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